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Growing evidence reveals universal hardening on various cosmic ray spectra, e.g. proton, positron, 
as well as antiproton fraction. Such universality may indicate they have a common origin. In this 
paper, we argue that these widespread excesses can be accounted for by a nearby supernova remnant 
surrounded by a giant molecular cloud. Secondary cosmic rays (p, et) are produced through the 
collisions between the primary cosmic ray nuclei from this supernova remnant and the molecular 
gas. Different from the background, which is produced by the ensemble of large amount of sources 
in the Milky Way, the local injected spectrum can be harder. The time-dependent transport of 
particles would make the propagated spectrum even harder. Under this scenario, the anomalies of 
both primary (p, e~) and secondary (e*, p/p) cosmic rays can be properly interpreted. We further 
show that the TeV to sub-PeV anisotropy of proton is consistent with the observations if the local 
source is relatively young and lying at the anti-Galactic center direction. 


PACS numbers: Valid PACS appear here 


I. INTRODUCTION 


It is widely accepted that cosmic rays (CRs) below the 
knee region originate from the galactic supernova rem- 
nant (SNR). In light of both primordial diffusive shock 
acceleration(I B] and steady-state transport|[3, Zi pro- 
cesses, the galactic CRs are expected to fall off as a fea- 
tureless power-law in a wide energy range from tens of 
GeV to ~ PeV. 

However, such a simple picture has been challenged by 
new observations in the past few years. The anomalies 
of CR electrons/positrons have been observed by quite a 
number of experiments, such as HEAT [5, e, AMS-01[7], 
ATIC{gh, PPB-BETS|9], PAMELA|10], Fermi-LAT|L1}, 
and most recently by AMS-02(12, . The overabun- 
dance of positrons requires extra primary sources, which 
involve either astrophysical objects, such as ulsars[14)- 
and the hadronic interactions in SNRs 2d 29 or 
more exotic origins like the dark matter self-annihilation 
or decay [23 H28}. For an extensive introduction of relevant 
models, one can refer to the reviews[2JH33} and references 
therein. 

Observations of proton, helium, and heavier nuclei 
also show remarkable hardenings at energies above a 
few hundred GeV/nucleon . The hadronic hard- 
ening may be ascribed to source model|4046], transport 
effect |47}152], or production mechanism|53, |54]. Most re- 
cently, AMS-02 collaboration released the measurements 
of antiproton-to-proton ratio, which shows a flat be- 
havior up to ~ 400 GeV [55]. This is not expected by 
the conventional production and propagation pattern of 
antiprotons|56, . Although sizable errorbars at high 
energies makes this conclusion not overwhelmingly yet, 
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it is worth to chew over the physical implication of a flat 
p/p fraction|58}{63}. 


In this paper, we aim at relating all these anomalies 
together in a unified framework. We propose that a local, 
fresh SNR, which was surrounded by a giant molecular 
cloud(MC), can explain all the above mentioned phenom- 
ena. In this scenario, both primary nuclei and electrons 
were accelerated by the shock wave in the local SNR, 
whereas the redundant positrons and antiprotons were 
yielded through the hadronic interactions between the 
primary nuclei and molecular gas. Compared to the back- 
ground component from the ensemble of large amount of 
Galactic SNRs, the injection spectra of this local SNR is 
harder. Similar picture has been suggested by [2], led, 
in which they focused only on the secondaries, namely 
positrons and antiprotons. Here we further extend the 
model to study its influence on primary CRs. Moreover, 
the consequence on Boron-to-Carbon (B/C) ratio is eval- 
uated in a realistic way, with reasonable assumption of 
the source nuclear abundance. This is especially impor- 
tant because the preliminary data published by AMS-02 
collaboration do not show significant deviation from the 
prediction of standard propagation modell56]. This, in 
turn, can impose strong constraints on the models which 
generate abundant secondary particles[65). Thirdly, we 
deduce both background and local SNR parameters co- 
herently without assumptions of the background param- 
eters a priori, which renders our treatment more self- 
consistent. 


The rest of paper is organized as follows: in section 
we briefly review the production and propagation of 
CRs in our model. The results and some discussions are 
presented in section [ILI Finally in section we give our 
conclusion. 
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II. COSMIC RAY TRANSPORT MODEL 


After being injected into the interstellar space, CRs 
diffuse within the Galactic magnetic halo by scattering 
off magnetic waves and MHD turbulence. The magnetic 
halo is usually approximated as a cylinder, with its radial 
boundary equal to the Galactic radius R = 20 kpc. Its 
half thickness L, which characterizes the vertical stretch 
of interstellar magnetic field, is fixed by CR data. Both 
CR sources and interstellar medium (ISM) are chiefly 
distributed in the Galactic disk, whose average width is 
roughly 200 pc, much less than halo’s thickness. The 
transport process of CRs in the magnetic halo is de- 
scribed by the following diffusion equation|3, 4, i.e. 


by a 3 Y 
Dr = Q(r,p) +V. (Drs VY E Vow) + Op PP S| 
o 
-5 [ee - 507 -vow]- 2-2. @) 


Here w(r,p,t) = dn/dp is the CR density per total par- 
ticle momentum p at position r. At the halo bound- 
ary, free escape condition is applied by default, namely 
U(R,z,p) = W(r,£L,p) = 0. The diffusion coefficient 
Dy» is supposed to be isotropic on large scales, which 
grows with particle rigidity R = pc/Ze, viz., 


R 


Dee = D9 (EY, (2) 


where 8 is the particle velocity in unit of the speed of 
light c. Both Do and 6 are specified as free parameters. 
Apart from diffusion, CR particles may still suffer from 
galactic convection V}, diffusive reacceleration Dpp, frag- 
mentation Tf, radioactive decay Tr, energy loss p, etc. 
In this work, we adopt the diffusive-reacceleration(DR) 
scenario |66]. The diffusive reacceleration originates from 
the random motion of magnetic field in ISM, which brings 
about the second-order Fermi acceleration during the 
transport. It is usually depicted as the diffusion in mo- 
mentum space, whose diffusion coefficient Dpp is related 
to spatial diffusion coefficient D,, and Alfvénic velocity 
va by 


Ap*v4, 
DppDrs = ————— A. 
pp 35(4 — 02)(4 —0) (3) 


In this work, the Galactic SNRs are separated into 
two classes, one is the local fresh SNR and the others 
are background ones. For the background SNRs, it is 
reasonable to assume that the spatial distribution of CRs 
from them arrives at steady state. Nevertheless for the 
local single SNR, the time-dependent transport of CRs 
after injection is requisite. In the next subsections, we 
discuss these two components separately. 


A. Background Supernova Remnants 


All of the SNRs other than the nearby one are labeled 
as background sources. Here for simplicity, the spatial 
distribution of background SNRs is assumed to be ax- 
isymmetric, that is, 


jea (E) wo [4 8 of 


with rg = 8.5 kpc, the distance from the solar system to 
the Galactic center. The parameters a and 8 adjusted to 
be compatible with the Fermi-LAT gamma-ray datal67]. 
The injection spectrum of both primary CR nuclei and 
electron is paramterized as a broken power-law, namely 


R\"% 
x br 
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Some other species, e.g. Li, Be, B, e+, p and radioac- 
tive elements, are hardly synthesized during the stellar 
nucleosynthesis. They can be brought forth from the 
fragmentation of parent nuclei throughout the transport. 
For the production of Li, Be and B, usually the so-called 
straight-ahead approximation is widely used, in which 
the kinetic energy per nucleon is conserved during the 
spallation process. The production rate is thus 


Qj = 


(5) 


XO (MUTH + mHetitHesj)vvi, (6) 
i=C,N,O 


where ny/He is the number density of hydrogen/helium 
in the ISM and oj4H/He-+,; is the total cross section of 
the corresponding hadronic interaction. 

Unlike above secondary CR nuclei, secondary e* and 
p have energy distribution. Therefore the source term of 
both et and p is the convolution of the energy spectra of 
primary nuclei ®;(£) and the relevant differential cross 
section doj4H/He+;/dE;, ie. 


Oe D [eve {rm SPP) 


dp; 
i=p,He Pj 


tnye tenali B | Gap), e) 
Pj 

Furthermore, antiprotons may still undergo non- 
annihilated inelastic scattering with ISM protons dur- 
ing propagation, in which antiprotons lose a significant 
amount of their kinetic energies. This is also known as 
the tertiary production. The numerical package GAL- 
PROH] is introduced to solve the transport equation 
to obtain the background CRs. The transport parame- 
ters are fine-tuned to fit the data together with the con- 
tribution from local source. 


1 http://galprop.stanford.edu/ 
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B. Nearby Young Supernova Remnant 


We assume a local(< kpc) supernova explosion oc- 
curred in a giant MC about 10° — 10° years ago. 
The charged particles were continually accelerated near 
the shock front with the expansion of supernova 
ejecta. The accelerated spectrum is represented by 
a power-law plus an exponential cut-off, ie. Q; = 
q;(R/Ro) * exp|-E/Ecut]. The normalization q; for 
each kind of element relies on the local chemical envi- 
ronment. In this work, we suppose that the element 
abundances in the MC are identical with the average of 
Galaxy. 

Besides, the CR nuclei generated by the local SNR also 
collide with the molecular gas around them and give birth 
to prolific daughter particles, like B, e+, p, and so forth. 
The yields of B and e*, p inside the MC are respectively 


Qj = 5 (NUCi+H—j + NHeCi+He—j VQ (Eteo (8) 


i=C,N,O 
and 
+00 d 
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where ny/He is the number density of hydrogen/helium 
in MC. In this work, we assume that it is 1000 times 
greater than the mean value of ISM. teol is the duration 
of collision. Q;(£) is the accelerated spectrum of primary 
nuclei inside local SNR. 

When the radius of shock front is comparable to or 
even larger than the size of MC, the whole MC was fi- 
nally fragmented subject to the expansion of ejecta. All 
the CRs break out of the MC and diffuse into the in- 
terstellar space. The time-dependent distribution of CR 
nuclei from a point source can be obtained by means 
of Green function technique. Since the local SNR only 
makes significant contribution above tens of GeV, a sim- 
plified transport equation can competent here with other 
terms, for example reacceleration, energy loss, etc., ne- 
glected. Then the transport equation for nuclei is 
= Q;(E)6(r—r')6(t—-t’)—2hd(z) P59; , 

(10) 
Here the transport equation is rewritten in terms of 
p = dn/dE with E the energy per nucleon. 2hd(z)I'; is 
the fragmentation term and TJ, = (nHO;+H+NHeT;}+He)U 
denotes the spallation rate of CR nuclei 7. Meanwhile ow- 
ing to that the interstellar medium is chiefly concentrated 
in the Galactic disk, which is much less than halo bound- 
ary, thin-disk approximation is applied. The correspond- 
ing analytical solution can be found in uA. As for the 
energetic electrons and positrons, the energy loss term, 
due to synchrotron in interstellar magnetic field and in- 
verse Compton scattering off CMB photons, is far more 


important here compared with fragmentation. Thus the 
transport equation for electron and positron becomes 


Fi _ D.. Ags + (Èp) = Q;(B)5(r — r')5(t — t') « 


(11) 
The energy loss rate E is approximately written as E = 
—bE? in Thomson limit. Here we use the analytical so- 
lution available in 


II. RESULTS 


In this paper, we apply our model to simultaneously 
fit the spectra of p/p, p and e*, as well as B/C and 
e7. The local SNR is supposed to lie at 0.1 kpc away 
from our solar system. In our propagation model, the 
essential transport parameters are Do, 6, va and L. The 


source parameters of the background CRs are Ap, vy, vy 


for proton, and A ,v§ ,v§ for electron. We still al- 
low for the fluctuation of positron background, which is 
denoted by a multiplier cet [66]. For the local SNR, the 
proton and electron parameters are respectively qj, aP, 
EP, and g§, a°. Moreover, we have another parameter 
specifying the collision duration of CRs within the MC, 
teol: For the heavier nuclei, like C, N, O, they share the 
same power-law index of injection as proton. In addition, 
to fit the low energy data, we have to take solar modula- 
tion into account, which is represented by the modulation 
potential o3]. The parameters of transport, background 
sources and local SNR are summarized in the separated 


Tables and 


A. Energy spectra 


Figure [1] illustrates our fittings assuming a compara- 
tively young local SNR, whose age is 10° year old. In 
the figure, the fluxes/ratios from the background(green 
dash), the local primary protons and electrons(blue dash- 
dot), the local secondary electrons and positrons(brown 
dash-dot) and total(black solid) are respectively repro- 
duced. Both background and total fluxes/ratios have 
been solar-modulated to account for the low energy spec- 
tra. The modulation potential for p(p), e+ and B/C are 
respectively ¢ = 660, 1300 and 330 MeV. 

For B/C ratio, it can be seen that both background 
and total can well fit the AMS-02 data. In this energy 
range, the contribution from local SNR is negligible, and 
hence the differences between two of them are tiny. The 
required transport parameters are Dp = 5.9 x 1078 cm? 
st, 6 = 0.343, va = 30 km s7! and L = 4.85 kpc, 
which are similar to other fittings However above 
600 GeV, the local component starts to dominate which 
renders the total B/C ratio hardening. This character 
can be verified by the future observations. Something 
similar also happens in the p/p ratio. From ~ 10 GeV, 
the total p/p ratio gradually deviates from the prediction 
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of conventional model and flattens up to 10? GeV, which 
agrees well with AMS-02 data. 

Since both positrons and anti-protons are mainly gen- 
erated by CR protons, the spectra of proton, p/p and 
e* collectively pin down the components of background 
and local proton. We find in order to fit these spectra all 
together, the power-index of local protons œ = 1.837 is 
harder than the background v2 = 2.49. The high energy 
cutoff Eeut of local proton is adjusted to 1.4 x 10° GeV 
for the sake of fitting high energy positron data. And to 
produce enough positrons at high energy, the duration of 
collision is 3.5 x 10? years when the density of local MC 
n = 1000 cm-?. 

For electron spectra, the local component comes from 
two origins, one is the shock-accelerated primary elec- 
trons within local SNR and the other is the products of 
the collision between CR nuclei and molecular gas. The 
normalization and injected power-law index of local pri- 
mary electrons are respectively g§ = 104° GeV~! and 
a = 2. Therefore the electron to proton ratio in local 
SNR is smaller than Kep = 0.01. And the injected spec- 
trum is also harder than background electron vg = 2.83. 
Figure[2]shows the similar results but hypothesizing an 
older local SNR(DR-B), whose age is fitted to be 3 x 10° 
years old. Both fitted transport and background param- 
eters are close to DR-A, but the injected power of local 
SNR has to be boosted to qo = 95 x 104° GeV~!, nearly 
seven times higher than DR-A. The total energy spectra 
of B/C as well as p/p resembles the model DR-A. Set- 
ting against Figure [1| the major differences come from 
the spectra of electron and positron, in which the total 
flux plunges dramatically at lower energy, around hun- 
dreds of GeV. Because of earlier releasing time of local 
CRs and energy loss during their propagation, only lower 
electrons are observed at present. This can be validated 
by the DAMPE experiment [69]. 


B. Anisotropy 


As exemplified above, both local SNRs with different 
ages can well explain the current data. In order to dis- 
tinguish the two cases, we further compute their accom- 
panied anisotropies of protons and electrons. The dipole 
anisotropy is defined as 


_ 3Dre |V] 
2 


ô (12) 


Owning to the spatial distribution of Galactic SNRs, 
there is inevitably a radial gradient of CR density, whose 
direction is toward Galactic edge. Therefore in the sce- 
nario of steady-state propagation, the anisotropy grows 
with diffusion coefficient, namely rises with the energy. 
However this is incompatible with current observations, 
e.g. proton(7Q], which do not manifest obvious energy- 
dependence between 1 TeV and 100 TeV. One of the in- 
terpretations is that there is a fresh nearby SNR(or more) 


lying at the anti-galactic center, whose CRs can offset the 
background streaming. 

In Figure[B] we show the anisotropies of protons respec- 
tively in DR-A(left) and DR-B(right). The blue dash-dot 
line is evaluated from the background and the black solid 
line is the total one while adding the local SNR. Appar- 
ently, the nearby source at the anti-galactic center direc- 
tion can effectively lower the anisotropy, but the extent of 
drop-off depends on the local SNR’s age. For the younger 
one, the total anisotropy can well conform with the mea- 
surements from 1 TeV to tens of TeV. But above 100 
TeV, due to the high energy cut-off of local proton flux, 
the anisotropy comes back to the case of steady-state. 

Due to the relatively close origin, the discreteness 
of the sources is more important for the anisotropy of 
electrons. However for the electrons, current measure- 
ments only give the upper bounds[7il FA. We com- 
pute the influence of anisotropies under both local source 
models, as shown in Figure [4] The blue dash-dot and 
black solid lines are respectively the anisotropies from 
the background assuming continuous distribution of the 
background sources and the total. Compared to DR-A, 
the local streaming in DR-B is comparable to the back- 
ground and the anisotropy at hundreds of GeV is tremen- 
dously suppressed. Nevertheless both anisotropies are 
still much less than the current constraints from 1-year 


Fermi-LAT|7i] and AMS-02[7). 


IV. CONCLUSION 


Currently, the excesses seem prevalent in both primary 
and secondary galactic CRs, according to the recent ob- 
servations. In this work, we address ourselves to this uni- 
versal anomaly in a unified viewpoint. Specifically, we 
envisage there was once a nearby supernova explosion, 
which occurred within the giant molecular cloud. The 
primary CRs are accelerated by the shock wave. And 
the secondaries, e+, p and B are massively produced by 
the interactions between primary CRs and molecular gas. 
All of them get free and enter into the interstellar space 
once the molecular cloud is fragmented. 

In order to fit the spectra of proton, p/p ratio and 
positron, the power index of local source has to be harder 
than the average of background SNRs. For the B/C ratio, 
the low energy spectrum is unaffected by the local source. 
However beyond 1 TeV, the total B/C ratio is well above 
the background. This is different from the spectrum of 
p/p ratio, whose transition occurs at lower energy. This 
stems from the differences of production cross section be- 
tween boron and antiproton. Future measurements ex- 
tending to higher energies can verify this notable feature. 

For the local SNR, we hypothesize two different ages, 
both of which can well recreate all the spectra. The main 
difference is at the high energy cut-off in the spectra of 
electron and positron. For the younger SNR, the cut-off 
occurs above 1 TeV, while for the older, the spectrum 
falls off well before 1 TeV, about hundreds of GeV. The 
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Fig. 1. The fittings to CR energy spectra, B/C(top left), p/p(top right), proton flux(middle left), positron flux(middle right) 
and electron flux(bottom left), by DR-A model. For B/C ratio, both AMS-02(red) and ACE(1998/01-1999/01) (blue) [za] data 
are used. The blue dash-dot is the flux from local SNR. The green dash and violet solid lines are respectively the background 
and total flux(or ratio) after solar modulation. The parameters of transport, background and local SNR. are respectively 
reproduced in table [I] [O] anda II] For B/C, p and e*, the parameter of solar modulation op/c = 330 MeV, ¢5/p = 660 MeV 
and @,4 = 1300 MeV. 
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Fig. 2. The fitting to CR energy spectra, B/C(top left), p/p(top right), proton flux(middle left), positron flux(middle right) and 
electron flux(bottom left), under DR-B model. For B/C ratio, both AMS-02(red) and ACE(1998/01-1999/01) (blue) ira) data 
are used. The blue dash-dot is the flux from local SNR. The green dash and black solid lines are respectively the background and 
total flux(or ratio) after solar modulation. The parameters of transport, background and local SNR are respectively reproduced 
in table[I] [O] and LI] For B/C, p and e*, the parameter of solar modulation op/c = 330 MeV, p = 600 MeV and ¢,+ = 1300 
MeV. 


model Do [10*8cm?/s] ô 
DR-A 5.9 
DR-B 5.5 


0.343 
0.383 


Ro [GV] va [km/s] L [kpc] 


4 30 4.85 
4 30 4.85 


Tab. I. Set of transport parameters. 


model v? RP. [GV] v} EP, [GeV] A5 v? 


DR-A 1.81 
DR-B 1.70 


12.88 2.49 100 
12.88 2.46 100 


Rf, [GV] v$ Ego [GeV] A8 


4.16 1.79 2.74 2.83 25 15.4 
4.05 1.70 2.74 2.88 25 


14.45 


Tab. II. Sets of parameters of background CR injection and solar modulation. The normalization flux of proton A$ is in units 
of 107° / (cm? - sr - s - MeV) and the electron’s A is 107'°/(cm? -sr - s - MeV). 


model aP 
DR-A 1.837 14 12.88 1.4e5 
DR-B 1.837 100 12.88 1.4e5 


q5 [10° GeV~"] Ro [GV] EP, [GeV] a® q6 [10° GeV~"] d [kpc] tage [10° yr] tco [10? yr] 


0.1 0.1 1 3.5 
1 0.1 3 3.5 


Tab. III. the parameters of local young SNR. a? and a° are the injected power-index for proton and electron. q and gj are 
respectively the source term of proton and electron at rigidity Ro = 12.88 GV. d is its distance to solar system and t the 
injection time of CR. teol is the interaction time between CRs and molecular gas. 


ongoing DAMPE experiment [69] could supply more pre- 
cise measurement at this energy range. 

For the local source model, one of the validation meth- 
ods is to compare its anisotropy of high energy CRs 
with observations. Under steady-state assumption, the 
anisotropy is proportional to the diffusion coefficient, 
which is in conflict with the accumulative data. If a local 
SNR is located at anti-Galactic center direction, the CR 
flux from it could effectively counteract the streaming 
from the Galactic center. We compute the anisotropies 
under our local SNR models. We find that when the lo- 
cal SNR is younger, the total anisotropy of proton is in a 
good agreement with the observations from 1 TeV to tens 
of TeV. In the case of electron, the local streaming can- 
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